Hysteretic memory effects in disordered magnets 
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We study the return point as well as the complementary point memory effect numerically with 
paradigmatic models for random magnets and show that already simple systems with Ising spin 
symmetry can reproduce the experimental results of Pierce et al. where both memory effects be- 
come more pronounced for increasing disorder and return point memory is always better than 
complementary point memory. 
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Hysteresis is ubiquitous in nature: it occurs in inter- 
acting systems as a collective phenomenon, in most disor- 
dered systems, across first-order transitions, in magnetic 
systems, and in depinning phenomena, among others. It 
is also crucial for industrial applications: most notably 
hysteresis lies at the very foundation of the magnetic 
recording industry. Hysteretic systems are employed as 
recording media because they retain their state for a long 
period after a write operation: they exhibit memory. 
Many aspects of this memory effect have been studied: 
the long-time decay of information through the super- 
paramagnetic decayi or how to optimize the signal-to- 
noise ratio by magnetic modeling, among others. 

Because of the central importance of memory in record- 
ing media, other types of memory effects have been stud- 
ied as well. Pierce et al. have investigated Co/Pt multi- 
layer films by repeated cycling over the hysteresis loop. 2 
The macroscopic magnetization returns to the same value 
cycle after cycle, a signature of a further type of mem- 
ory. It is a profound question whether this memory comes 
about by the system returning to the same microscopic 
configuration: an effect dubbed "return point memory'' 
(RPM), or whether it is exhibited only on the macro- 
scopic level. Up to the groundbreaking work of Pierce 
et at, it was not possible to address this question ex- 
perimentally. By using an x-ray magnetic speckle micro- 
scopic technique which can indirectly resolve microscopic 
domain patterns, Pierce et al. found the presence of RPM 
for strong disorder, but no RPM for weak disorder. Fur- 
thermore, Pierce et al. studied a second type of memory: 
whether the system develops the reverse of the original 
configuration when the hysteresis sweep reaches the op- 
posite magnetic field, called "complementary point mem- 
ory" (CPM). They found that disorder influences CPM 
in the same way as RPM. In addition, they report that 
CPM appears to be smaller than RPM. However, the 
RPM-CPM difference does not exceed 10% and it is not 
clear whether the difference could have been caused by 
instrumental bias. 

The experimental work of Pierce et al. spawned two 
theoretical studies to date. Deutsch and Mai have per- 
formed micromagnetic simulations using the Landau Lif- 
shitz Gilbert (LLG) equations^ reproducing the above 



memory effects. However, only few disorder realizations 
were considered, leading to large error bars. Further- 
more, the difference between RPM and CPM was not 
much bigger than the error bars of the simulation. Fi- 
nally, they suggest that the rotation of spins is primar- 
ily responsible for the RPM-CPM difference, i.e., scalar 
spins cannot cause the effect (in the absence of random 
fields). But since their model Hamiltonian contained nu- 
merous terms with several parameters, further clarifica- 
tion may be useful to determine which of these terms is 
the primary cause of the memory effects. Jagla also used 
the LLG approach, but for constrained soft scalar spins, 
subject to random fields. 4 He reported seeing the memory 
effect and the RPM-CPM difference as well. However, 
no quantitative measure of the memory was given and 
the simulations were performed only at T = 0. 

Therefore, the following challenges still remain to be 
understood: (i) What is the primary cause of the mem- 
ory effect, or, equivalently, what is the minimal model 
that exhibits the memory effect? (ii) Does the memory 
effect persist at finite temperatures, since the thermal 
fluctuations have the potential to destroy microscopic 
correlations? (iii) What is the disorder dependence of 
the memory effects? (iv) Does the RPM-CPM differ- 
ence convincingly exceed the error bars? To address 
these challenges, we have studied RPM and CPM in min- 
imal, paradigmatic disordered spin models: spin glasses 
and random-field models. Equilibrium and close to equi- 
librium dynamical properties of spin glasses have been 
studied extensively over the years^&i However, far-from- 
equilibrium properties, such as their hysteresis, are much 
less understood. Recent studies characterized the hys- 
teresis loop of the Edwards- Anderson spin glass (SG), 8 
the random field (RF) Ising model, 9 and the Shcrrington- 
Kirkpatrick models While the RF model has been stud- 
ied before, a comprehensive study with different models 
is still lacking. In this paper we address the above four 
challenges by studying disordered Ising-type models with 
and without frustration. 

We first study the nearest-neighbor Edwards- Anderson 
Ising spin glass 5 with the Hamiltonian 
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in which Ising spins Si = ±1 lie on the sites of a square 
lattice of size N = L 2 with periodic boundary conditions. 
The measured quantities show essentially no size depen- 
dence past L — 20. The interactions Jij are Gaussian dis- 
tributed with zero mean and standard deviation aj . The 
simulations are performed by first saturating the system 
by applying a large external field H and then reducing 
H in small steps to reverse the magnetization. At zero 
temperature we use standard Glauber dynamics^, where 
randomly chosen unstable spins (pointing against their 
local field hi = Yl{j) Jij^j +-^0 are flipped until all spins 
are stable for each field step. At finite temperatures we 
perform Monte Carlo simulations until the average mag- 
netization is independent of the equilibration time of the 
Monte Carlo simulation for each field step. All results 
are averaged over 500 disorder realizations. 

Pierce et al. captured the RPM and CPM in terms of 
overlaps of the spin configurations at different points of 
the hysteresis loop. Accordingly, we capture CPM and 
RPM with q(H*), the overlap of the spin configuration 

sf^ at a field H* with the configuration S^ at a field 
with the same magnitude \H*\ after n = l/2 (CPM) and 
1 (RPM) cycles around the hysteresis loop, respectively: 




FIG. 1: (a) In CPM spin configurations on opposite branches 
I and II [n — 1/2 in Eq. °f the hysteresis loop (gray dots) 
are compared, whereas in RPM a configuration on branch 
I [{£} } in Eq. is compared to itself after n £E N cycles 
around the loop (black dot), (b) Configurations for \H*\ = 1.2 
for the two-dimensional Edwards- Anderson SC. Light pixels 
correspond to down spins, dark pixels correspond to up spins. 
Top: initial configuration at H* — —1.2 (branch I). Center: 
configuration after a half cycle at H — —H* (CPM, branch 
II). Bottom: configuration at H — H* after one cycle around 
the hysteresis loop (RPM). The black pixels represent differ- 
ences between the initial and final configurations. 



q(H*) = [ -ZJ_Y,S? ) {H*}st ) {(-l) 2n H*} . (2) 

i=l 

The uniqueness of CPM and RPM is tested by q'(H, H*), 
the overlap between the spin configuration at H* 

with a configurations s!- at a field H after n = 1/2 
(CPM) and 1 (RPM) cycles around the hysteresis loop: 

/ \2n ^ 

q'(H,H*) =SZJ_J2S^[H*]S^{H] . (3) 

i=l 

To develop a physical picture, Fig. illustrates configu- 
rations relevant for the CPM and RPM memory effects. 
The SG displays a large degree of CPM/RPM. 

Figure |21 (left panel) shows q(H*) as a function of H* 
for various temperatures. The RPM and CPM curves 
are indistinguishably close. The strong CPM can be at- 
tributed to the spin-reversal symmetry of the system: 
upon reversing all spins Sj and the magnetic field H, 
the Hamiltonian transforms into itself. Note that RPM 
is not perfect at T = because of the nature of our spin 
updating. If unstable spins were sorted and not picked 
at random, RPM would be perfect at T = (this is not 
the case if the ground state of the system is degenerate). 
The observation of a robust RPM and CPM memory is 
an answer to challenge (i), establishing the SG as a pos- 
sible minimal model displaying memory effects just as in 
the experiments of Pierce et al. with no adjustable pa- 
rameters. Also noteworthy is that the CPM and RPM 
are smallest around the (temperature-dependent) coer- 
cive field: as the number of equivalent spin configurations 
is the largest in that field region, the reversal process can 
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FIG. 2: (Color online) Overlaps for the 2D Edwards- Anderson 
SG. Left: Overlap q [Eq. @] for different values of H* at 
different temperatures T. Right: Overlap q [Eq. Jjjj] as a 
function of the applied field H for H* — —1.2 (coercive field 
H c — —0.98). The memory becomes more pronounced for 
T — > 0. Data for aj = 1. The dotted line represents the 
zero-temperature major loop. 

evolve along many different paths. We also address the 
temperature dependence of the overlap, as raised at (ii). 
FigureHshows that the T = CPM and RPM survive to 
finite T, even though it would be natural for the thermal 
fluctuations to wash out the memory at the microscopic 
level and convert it to a macroscopic memory only. In 
addition, the memory decreases with increasing temper- 
ature. 

Figure |21 (right panel) shows q'(H,H* = —1.2) for 
various temperatures. The data show the uniqueness 
of CPM and RPM: the overlap function strongly peaks 
at H = \H*\. Thus, the memory is not the result of 
a gradual slow buildup: there are large-scale spin rear- 
rangements during the field sweep in frustrated systems, 
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FIG. 3: (Color online) Left: overlap q(H*) for RPM and CPM 
(2D SG) as a function of H* for different disorder strengths 
a j. Right: minimum value g m i n (at coercivity) of q(H*) as 
a function of disorder for RPM and CPM. Both sets of data 
show that the RPM and CPM increase with increasing disor- 
der. 

recreating the initial spin configuration only in the close 
vicinity of H = \H*\. As above, CPM and RPM decrease 
with increasing temperature. The nonzero plateaus are 
caused by the measuring field H* = — 1.2 being different 
from the coercive field, where M(H*) ^ 0. 

Next, we address (iii), i.e., we explore the disorder de- 
pendence of both memory effects, since this aspect was 
an important element of the experiments of Pierce et 
al. Figure 13 (left panel) shows q(H*) as a function of 
the disorder aj at T = 0.2. The right panel of Fig. 
shows q(H*) at coercivity as a function of oj. RPM and 
CPM are indistinguishable and visibly increase with in- 
creasing disorder, in good agreement with Ref. 0. The 
physical reason for this correlation is that for weak dis- 
order, the energy landscape includes many comparable, 
shallow valleys, without a single optimal path. There- 
fore, during subsequent field sweeps, the system evolves 
along different paths, thus reducing the memory. In con- 
trast, for stronger disorder, the energy landscape devel- 
ops a few preferable valleys and the system evolves along 
these optimal valleys during subsequent cycling around 
the hysteresis loop. [The dips in q(H*) become broader 
with increasing disorder because the entire hysteresis 
loop broadens.] Finally, in relation to (iv) it is noted that 
in the SG RPM and CPM are indistinguishably close. 

Next, we explore RPM and CPM, as well as the same 
four challenges in the 2D random field Ising model, 

where J = 1 is a ferromagnetic coupling and the random 
fields hi are chosen from a Gaussian distribution with 
zero mean and standard deviation ah- The main differ- 
ences between the RF model and the Edwards- Anderson 
SG are that the RF model does not have frustration and 
does not have spin-reversal symmetry. 

Figure H (left panel) shows both q(H*) for the RPM 
(dashed lines) and the CPM (solid lines) as a function 
of the field H* for various temperatures. Concerning 
(i), this model also exhibits clear RPM, whereas CPM 
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FIG. 4: (Color online) Left panel: Overlap q for RPM (dashed 
lines) and CPM (solid lines) for different temperatures for the 
RF model. The RPM is perfect at T = 0. For all tempera- 
tures the RPM is better than the CPM. In particular, for all 
temperatures the CPM is close to zero at coercivity suggest- 
ing that CPM is rather poor in this model. Right panel: Field 
dependence of q{H*) for RPM for different disorder strengths 
ah (T = 0.20). The RPM becomes better with increasing 
disorder at finite temperatures. The inset shows the disorder 
dependence of q(H*) for the RPM at coercivity. 



is rather poor. Further, (ii) RPM is not washed out by 
thermal fluctuations at once, but is gradually weakened 
with increasing T. 

Regarding (iv), the RF model deviates from the SG 
results and correlates with the experiments: in the RF 
model the RPM and CPM are different. The RPM is big- 
ger than the CPM for all temperatures. This effect is due 
to the fact that the RF model does not have spin-reversal 
symmetry and therefore the spin configurations on the 
ascending branch do not correlate closely with the con- 
figurations on the descending branch. For intermediate- 
to-large values of the disorder the CPM is negligible and 
in the proximity of the coercive field the CPM correla- 
tion is in fact negative. In contrast, the RPM is large in 
the RF Ising model. In particular, at T = the RPM is 
perfect due to the "no-crossing property" of the RF Ising 
model. 11 Figure 01 (right panel) shows (iii) q(H*) for the 
RPM for different disorder strengths a h at T = 0.2. RPM 
in the RF model also increases with increasing disorder, 
also illustrated in the inset: The overlap q(H*) at coer- 
civity for RPM increases with increasing disorder. As for 
the SG, the memory effects increase due to the valleys 
in the energy landscape becoming more pronounced with 
increasing disorder. 

In the RF model the CPM is always close to zero, 
whereas the RPM increases with increasing disorder (and 
is perfect at T = 0). Correspondingly, the RPM-CPM 
difference is large over much of the parameter space. 
These findings do not agree completely with the experi- 
mental results of Pierce et al. On the other hand the SG 
lacks the RPM-CPM asymmetry completely. Therefore 
we have explored whether a combination of the SG and 
the RF model might yield results comparable to the ex- 
periments: increasing memory with increasing disorder, 
as well as RPM being better than CPM. Thus we apply 
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FIG. 5: (Color online) Left panel: Field dependence of q(H*) 
for RPM (solid lines) and CPM (dashed lines) for different dis- 
order strengths aj (T = 0) for the SG+RF model [Eq. ©]. 
The RPM is better than the CPM and both increase with in- 
creasing disorder. Right panel: Disorder dependence of q(H") 
of the RPM and the CPM at coercivity. Both memory effects 
become better with increasing disorder, and RPM is larger 
than CPM, in good qualitative agreement with the experi- 
mental results of Pierce et al. 

diluted random fields to the SG model: 

Hsg+RF = -J2 J » S * S 3 ~ E [ ff + W • ( 5 ) 

The random bonds Jy are chosen from a Gaussian distri- 
bution with zero mean and standard deviation aj. The 
random fields are randomly attached to only 5% of all 
sites, chosen from a Gaussian distribution with zero mean 
and standard deviation unity. 

Figure [S] illustrates the results at T = which persist 
well to finite T. The left panel of Fig. [3 shows q(H*) for 
both the RPM and the CPM as a function of H* for dif- 
ferent disorder strengths a,j in the SG+RF model. The 
right panel of Fig.[S]shows the overlap q at coercivity for 
the RPM and the CPM. Both memory effects increase 
with increasing disorder and the RPM is larger than the 
CPM, in qualitative agreement with experimental results. 
These findings establish that the SG+RF model qualita- 
tively reproduces all aspects of the Pierce et al results. 

In summary, we have simulated paradigmatic disor- 
dered spin models: the Edwards- Anderson spin glass, the 
random-field Ising model, and a spin glass with diluted 
random fields. We have found that (i) all three models 



exhibit return point and complementary point memory. 

(ii) Both memory effects persist to finite temperatures. 

(iii) Both memory effects increase with increasing disor- 
der, (iv) In the spin glass the RPM is always identical 
to the CPM because of the spin-reversal symmetry. In 
the RF Ising model the CPM is always close to zero be- 
cause of the lack of spin-reversal symmetry. Finally, a 
spin glass where spin-reversal symmetry is broken with 
diluted random fields shows a moderate RPM-CPM dif- 
ference, illustrating the sensitivity of the memory effects 
to the details of the models. 

In relation to the experiments of Pierce et al. , one re- 
calls that the films of Pierce et al. have strong out-of 
plane anisotropy, restricting the orientation of the spins 
to being essentially perpendicular to the film. Thus, de- 
scribing those films in terms of Ising Hamiltonians might 
be a reasonable approximation. Furthermore, since the 
dipolar interactions in these perpendicular films are anti- 
ferromagnetic, they introduce extensive frustration into 
the system, the key ingredient of the spin glass model. 
Finally, spins frozen in by shape anisotropics of a locally 
deformed environment, or by unusually large crystal field 
anisotropies, or possibly by frozen-in reversed bubbles, 
as reported by Davies et al. in Co/Pt multilayer films, 12 
could all be the origin of random fields at a few per- 
cent of the sites. These considerations make it conceiv- 
able that all ingredients of the SG+RF model may be 
present in the Pierce et al. films. In addition, the SG+RF 
model agrees qualitatively with the experiments: RPM 
and CPM increase with increasing disorder and RPM is 
always more pronounced than CPM. 

These results, together with the work by Deutsch and 
Mai, and that of Jagla show that a RPM-CPM asym- 
metry can only be obtained when the system's symme- 
try is broken. While Deutsch and Mai as well as Jagla 
break the symmetry dynamically using vector models 
with damped LLG dynamics, in this work we present 
an alternate avenue that breaks the symmetry statically 
using random fields. Thus in this work we present a 
plausible alternative which reproduces all aspects of the 
experiments within a minimal framework. 
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